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The Laser Interferometer Space Antenna (LISA) is expected to provide the largest observational 
sample of binary systems of faint sub-solar mass compact objects, in particular white-dwarfs, whose 
radiation is monochromatic over most of the LISA observational window. Current astrophysical 
estimates suggest that the instrument will be able to resolve ~ 10* such systems, with a large 
fraction of them at frequencies <!, 3 mHz, where the wavelength of gravitational waves becomes 
comparable to or shorter than the LISA arm-length. This affects the structure of the so-called 
LISA transfer function which cannot be treated as constant in this frequency range: it introduces 
characteristic phase and amplitude modulations that depend on the source location in the sky and 
the emission frequency. Here we investigate the effect of the LISA transfer function on detection 
and parameter estimation for monochromatic sources. For signal detection we show that filters 
constructed by approximating the transfer function as a constant (long wavelength approximation) 
introduce a negligible loss of signal-to-noise ratio - the fitting factor always exceeds 0.97 - for 
/ < 10 mHz, therefore in a frequency range where one would actually expect the approximation 
to fail. For parameter estimation, we conclude that in the range 3 mHz ^ / ^ 30 mHz the errors 
^H ' associated with parameter measurements differ from ~ 5% up to a factor ~ 10 (depending on the 

3 \ actual source parameters and emission frequency) with respect to those computed using the long 

wavelength approximation. 
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^v ' The Laser Interferometer Space Antenna (LISA) is a space borne laser interferometer of arm-length 5 million km 

(T^ \ for the observation of gravitational waves (GWs) in the frequency window 10^'* Hz - 0.1 Hz il|. Amongst the great 
^^ ■ variety of sources that the instrument will be monitoring jaL binary systems of faint sub-solar mass compact objects, 
^«0 ■ in particular white-dwarfs, will be considerably abundant H,|l, |M|y{. In fact LISA is expected to provide the largest 
^— ^ . observational sample of these faint stars; current estimates suggest that ~ 10^ systems will be resolvable |S], including 
about a dozen of already known galactic binary systems (the so-called "verification sources") easily detectable during 
the first few weeks of operation. Considerable attention has been recently devoted to (sub-) solar mass binary systems, 

^ ' as they will allow us to investigate the evolutionary history of degenerate stars, their formation rate, distribution in 
i^[ the galaxy and mass transfer [IjHiSH- Key issues in preparation of LISA are the design of appropriate data analysis 

^ • schemes to extract effectively and efficiently signals from noise and the investigation of the astronomical information 

Pv^] that can be gathered from the LISA data set, possibly followed up by observations with optical telescopes [lOl[lll |. 

^ . Sub-solar mass binary systems are expected to produce a moderate-to-large signal-to-noise ratio in the LISA data 

set and their radiation is monochromatic over most of the instrument sensitivity window, i. e. the intrinsic frequency 
drift during the observation time T w 1 yr is smaller than the frequency resolution bin of width A/ = l/T (this of 

C^ ■ course ignores the spreading of power in adjacent frequency bins induced by the LISA orbital motion around the Sun 
and the change of orientation of the detector with respect a putative source). In the context of data analysis, both 
signal detection and parameter estimation, the response of LISA to gravitational waves plays a vital role because it 
introduces features, such as amplitude and phase modulations, that need to be properly accounted for during signal 
processing. For LISA, which observes gravitational radiation from binary systems over a frequency band where the 
signal wavelength A can be either longer or shorter than the interferometer arm length, L = 5 x 10^ km, it is well 
known that the detector response changes dramatically at / ~ /* = 1/(2 ttL) ~ 9.6 mHz, which corresponds to the 
inverse of the round-trip light-time along the LISA arms [l^- If / <C /* (long- wavelength regime), the instrument 
transfer function can be effectively regarded as constant while for / ^ /* the transfer function depends on the source's 
instantaneous emission frequency and location in the sky, which in turn introduces time and frequency dependent 
amplitude and phase modulations at the detector output. The LISA transfer function, then, starts to substantially 
change the structure of the LISA Michelson observable right at the heart of the instrument's sensitivity window / ~ 3 
niHz 13J , and yet the implications for data analysis have been so far largely ignored. In the context of signal detection, 
no study has been devoted to the identification of the frequency range over which it is indeed safe (with respect to 
losses of signal-to-noise ratio) to approximate the transfer function as a constant in generating signal templates (which 
would clearly reduce the complexity of the problem at hand). As far as parameter estimation is concerned, a number 
of studies have been carried out so far for monochromatic sources aimed at computing lower bounds to the the errors 
associated with parameter measurements [sBISillS- However, in all of them the LISA transfer function has been 
treated as a constant. Such a simplifying assumption is likely to affect the evaluation of the actual accuracy with 
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which source parameters can be measured, and has already been pointed out in the context of observations of massive 
black hole binary systems p!5| . 

The goal of this paper is to investigate the effect of the LISA transfer function on signal detection and parameter 
estimation for monochromatic sources. In particular (i) we identify the frequency window over which signal templates 
can be safely constructed by approximating the transfer function as a constant (i.e. the signal-to-noisc ratio is very 
marginally affected) and (ii) explore the implication of the frequency dependent transfer function on the errors of 
parameter measurements. For signal detection, the main outcome of our work is that one can safely work in the long- 
wavelength approximation well into the regime where one would actually expect the approximation to fail: effectively 
the use of templates computed using the long-wavelength approximation does not introduce any significant loss of 
signal-to-noise ratio up to / « 10 mHz. For parameter estimation we conclude that at / ~ 5 mHz the errors are 
already starting to depart, by ss 5% — 10%, from those computed using the long wavelength approximation. Such 
discrepancy becomes more pronounced, on average, as the signal emission frequency increases, and in the frequency 
range 10 mHz < / < 30 mHz the errors are considerably smaller (up to a factor ~ 10 for the 3-armed LISA) than the 
ones previously reported in the literature. 

The paper is organised as follows: in Section E] we review the signal detected at the output of the LISA Michel- 
son interferometer in the so-called long-wavelength approximation (constant transfer function) and rigid adiabatic 
approximation (frequency dependent transfer function); Sections IIIII and IIVI contain the key results of the paper: in 
Section IIIII we estimate the degradation of signal to noise ratio (as a function of frequency) introduced by searching 
for signals using templates that approximate the transfer function as a constant; in Section llVI we show the effects 
of the LISA transfer function on parameter estimation by computing the inverse of the Fisher information matrix; 
in both cases we perform extensive Monte Carlo simulations in order to sample a wide parameter space; Section IVI 
contains our conclusions and pointers to future work. 

II. THE SIGNAL MEASURED AT THE LISA OUTPUT 

LISA consists of a constellation of three drag-free spacecraft placed at the vertices of an ideal equilateral triangle 
with sides L ~ 5 x 10^ km to form a three-arm interferometer, with a 60° angle between two adjacent laser beams. 
The barycentre of the instrument follows an almost circular heliocentric orbit (the eccentricity is < 0.01), 20° behind 
the Earth; the detector plane is tilted by 60° with respect to the Ecliptic and the instrument counter-rotates around 
the normal to the detector plane with the same period 1 yr (we refer the reader to [J, 0, Il3| and references therein for 
more details) . Due to the finite arm length, the round-trip light-travel time between two vertices of the constellation 
is finite and this introduces a characteristic frequency defined as: 

/* = ttV - 9-6 X 10"' Hz . (1) 

z ttL 

The LISA "Michelson output" is synthesised in software using a technique known as Time Delay Interferometry 
(TDI) by combining the Doppler readouts registered at several points in the LISA constellation 16] . The eccentricity 
of the orbit, the motion of the spacecraft and the finite arm length of LISA affect, in a non trivial way, the detector 
output generated by metric perturbations induced by impinging GWs (cf Eq. (ll)-(39) of [i3)- However, two fairly 
simple approximations to the exact LISA output have been derived: (i) the long wavelength approximation and (ii) 
the rigid adiabatic approximation. In both cases the eccentricity of the orbit and time dependency of the arms are 
ignored. In the long wavelength approximation the actual size of LISA's arms with respect to A (the wavelength of 
the impinging radiation) is also neglected, implying that the instrument transfer function is assumed to be constant. 
This is a very good representation of the LISA output for / <C /*. In the rigid adiabatic approximation the frequency 
and source position dependence of the transfer function are fully included; the LISA Michelson output that is derived 
under these assumptions is an excellent representation of the exact one up to / « 0.5 Hz, in the sense that the overlap 
between the exact and approximated output is always > 0.97 [13. In the remainder of the paper we will therefore 
consider the rigid adiabatic approximation as a faithful representation of the signal extracted at the LISA output, 
as we are considering binary systems whose radiation is at / < 0.5 Hz. We refer the reader to [l|, Q, U^ [13 and 
references therein for more details. Here we briefly review, mainly to establish notation, the expressions of the signal 
recorded at the output of the LISA Michelson interferometer using these two different approximations. 

Consider a generic gravitational wave source whose position in the sky, with respect to an observer on LISA, 
is identified by the unit vector N. Gravitational waves travel in the — N direction, and are described by the two 
independent polarisations /i+ and hx [13 ■ The metric perturbation hat at the detector can be decomposed as the 
sum of the two independent polarisation states according to 

Kbit)^h+it)e+,{t) + h-{t)e^,{t), (a, 5 =1,2, 3), (2) 



where e|^j and e^^ are the wave's polarisation tensors. They can be expressed as a function of the source basis tensors, 
e^j and e^j,, and the wave polarisation angle ip as: 



^afc = cos 2^-6+ -sin 2^-6^6, 
<b = sin2V'e+ +cos2V'e„^,. 



(3a) 
(3b) 



The source basis tensors are constructed by considering two unit vectors to" and n° orthogonal to each other and to 
the wave propagation direction — N, in order to form a left-handed Cartesian tern, according to: 
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niamb - UaUb . 



(4a) 
(4b) 



If the signal is described by the two polarisation amplitudes A^{t) and Ax{t), respectively, and the gravitational 
phase (j)Gw{t), then Eq. |(5J) becomes 



where we have explicitly included the dependence on time. 



(5) 



A. The rigid adiabatic approximation 



In the so-called rigid adiabatic approximation the strain at the detector induced by the metric perturbation, Eqs. (0 
and ©, is [13 



h{t)^D''\t)hab[t) 
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where D°'^{t) is the time dependent detector response tensor. If we consider the detector formed by the arms Z° and l^. 

- \j {j — 1, 2, 3) are the unit vectors along each of the LISA arms; in a frame attached to the solar system barycentre 
they are time dependent functions as a consequence of LISA's motion - the detector response tensor reads 
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Tj is the instrument transfer function, whose behaviour as a function of time and frequency is shown in Figures Q 
and 121 In general, this is a time dependent function which depends on the source location and the instantaneous 
frequency of the wave; it carries key information about the interferometer arm length through the characteristic 
frequency /*, Eq. (^. It is straightforward to check, cf Figures Q] and |21 that for ///* ^ 1, Tj tends to unity; as 
the frequency of GWs increases it develops increasingly stronger oscillations. If we ignore the contribution given by 
the geometry of the detector, the transfer function introduces a phase shift ^ 1 rad for / ^ 5 mHz (the LISA peak 
sensitivity is at / ~ 3 mHz). The overlap of the output © with the exact expression of the LISA readout is greater 
than 0.97 up to w 0.5 Hz |l3j . so that we can safely replace the full expression with the rigid adiabatic approximation 
in the frequency range of interest for this paper (0.1 mHz - 0.1 Hz). 

The expression for the detector response tensor, Eq. {Tj), can be recast in the form 
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FIG. 1: The time evolution of the LISA transfer function for selected frequencies. The plot shows the transfer function, cf 
Eq. (0 , as a function of time (for one LISA period, corresponding to 1 yr) for one of the arms and a randomly chosen source 
location. The values of the transfer function for four selected frequencies (see labels) are shown: / — 10"^ Hz, 10~^ Hz, 10~^ Hz 
and 10~ Hz (solid line: real part; dashed line: imaginary part). Notice that even by / = 10~ Hz the transfer function is no 
longer constant over the LISA orbital period. 
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and 



^1 = ^(3-^.^^^)' 
6 = ^(3-^fc^c), 
Is is also useful to introduce the generalised antenna beam patterns Fn and Fn defined as 
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FIG. 2: The frequency dependence of the LISA transfer function. The plot shows the LISA transfer function, cf Eq. Q, as 
a function of frequency for two (see labels) selected relative positions and orientations of LISA with respect to a source (solid 
line: real part; dotted-dashed line: imaginary part). Notice that at a few mHz the transfer function is already significantly 
diff'erent from the constant value apparent at / ^ /*. 



Fn and F„ are time dependent because of the change of orientation of LISA during the typical observation time 
~1 yr. It is straightforward to verify that at low frequencies 



F+ = J2fM if«f^). 



n=l 
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where F^ and F^ are the usual antenna beam patterns cf. Eqs. H2Ua|) and (|20b|l . 

Using Eqs |5)), ©, Hl()a(l - (|10d|l and (|lla|l - ljlld|l . the (real) strain at the detector becomes: 



(13a) 
(13b) 



h{t) = E(Fi+)(i)A+(i)cos[0GVv(t)-e„(t) + ^D(t)] 

n=l 

+Fi''\t) A^{t)sm[^Gw{t) ~ Ut) + ^B{t)]) , (14) 

where we have included the phase Doppler shift induced in the signal at the LISA output by the motion of the LISA 
barycentre with respect to a source, see Eq. (|A13p . Using double angle formulae one can express Eq. (|14|l as the 
superposition of 4 harmonics with different (time dependent) polarisation amplitudes and phases 
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and 



ifnit) = arctan 



(18) 



Fi+\t)A+{t) 
Explicit expressions for Bn{t), Xnit) and (pnit) are given in the Appendix. 

B. The long wavelength approximation 

In the low frequency region of the LISA sensitivity band, i.e. f / f* <^ 1, the transfer function can be approximated 
as a constant |13 , cf Figures ^ and |21 Tj is therefore independent of the signal frequency and source position in the 
sky. The strain at the detector output, Eq. ©, takes the usual form 

hL{t) - h+{t)F+it) + /ix (t)Fx (t) . (19) 

This expression is known as the long wavelength approximation to the detector output. In the previous expression F^ 
and i^x are the antenna beam patterns. If one takes the observable constructed using the two arms identified by Ij 
and Ifc, they read 

F+ = l{ip';-ltll)et,, (20a) 

Fx = l{lili--l2l2K,- (20b) 

As is the case for Fn and Fn' (n = 1, ..., 4), F^ and F^ are time dependent because of LlSA's change in orientation 
with respect to source location during the observation period. 
The detector output ifT^ can be cast in the form 

hL{t) = B{t) cosxit) (21) 

where 

B{t) =. [(F(+)(t)A+(t))2 + (F(x)(i)Ax(t))']'/', (22) 

X{t) = ^Gw{t) + ^{t)+(pu{t), (23) 



and 

Lp(t) = arctan 



(24) 



J^(x)(t)Ax(t) 
Fi+)it)A+{t) 

Explicit expressions for B{t), xi^) a-^d (f){t) are given in the Appendix. 

III. SIGNAL DETECTION 

In this section we investigate the impact of the instrument transfer function on detection. Due to the essentially 
perfect knowledge of signal waveforms for this class of sources, we assume that the signal processing scheme will be 
based on a coherent approach, where the data are correlated with a discrete bank of templates to extract the signal 
from the noise. The goal of this section is to identify the frequency band over which the use of the long-wavelength 
approximation Hl for signal templates (i.e. constant transfer function), cf Eq. H21|l. does not affect the signal-to-noise 
ratio (SNR) at which the actual radiation embedded in the detector noise can be detected. At the LISA output, in 
fact, gravitational waves are modulated by the complex structure of the instrument transfer function and are actually 
represented by Eq. ((TK|l . It has been recently suggested [iJl that the long wavelength approximation, ^.^(i), is already 
not a good approximation of the exact detector output at / « 3 mHz; here we show that for the purpose of detecting 
monochromatic signals h^^t) is actually a completely satisfactory approximation at least up to / ~ 10 niHz, in the 
sense that it returns a fitting factor which is always greater than 0.97. We start by briefly reviewing the key concepts 



and formulae of signal detection through matched filtering - we refer the reader to [3, ITa . Il9l | and references therein 
for more details ~ and then present the results of our analysis. 

The signal s{t) registered at the detector output is a superposition of noise n{t) and gravitational waves h{t; A), 

s{t)=h{t;\)+n{t), (25) 

where A represents the vector of the unknown parameters that characterise the waveform. In the case of monochro- 
matic sources h{t) is fully described by seven independent parameters: signal amplitude Aqw, frequency /o, ar- 
bitrary initial phase at the beginning the observation (f>Q, and four angular parameters related to the position 
of the source in the sky, 6]\[ and (f)]\r, and the orientation of the orbital plane, 9^ and 0l (we refer the reader 
to the Appendix where the explicit dependence of h(t) on A is presented). The parameter vector is therefore 
A = {AGw,4'o^foi(^N,4'N,0L,(j>L} = {^GW>o,/o,^}, where B = {6N,(f>N,0LAL}- Our notation reflects the fact 
that j4gw is not a search parameter, as it determines simply the signal-to-noise ratio at which detection is made, /o 
can be easily searched over using the Fast Fourier Transform of the data stream and 0o is an extrinsic parameter 
that can be trivially maximised. We assume the noise to be stationary and Gaussian (although we do not expect this 
condition to be fully met by the actual data), characterised by a noise spectral density <S'„(/). In the geometrical 
approach to signal processing the signal hit) represents a vector in the signal manifold and the signal parameters A 
are the coordinates on this manifold. One can introduce the following inner product between two signals v and w |7|: 

<«-/:: ^if^^^ 

= — / v*{t)w{t)dt, (26) 

where the second equality follows from Parseval's theorem and 5*0 is the (essentially constant) noise spectral density 
at /o. 

According to the definition 126|) . the optimal SNR at which h can be detected is 

S^ _ ih\h) _tu\uM/2 



.T ; r.., M -(/i|M ^ • (27) 

NApt rms[(/i|n)] 

If q{t] A) is the family of templates used to search for the class of signals h{t\ A) - notice that A, the template 
parameter vector, is not necessarily the same (including its dimensions) as A, and q does not necessarily belong to the 
same manifold as h - the adequacy of the template family q to search for the signal family h is given by the so-called 
fitting factor (FF), defined as 



quae 

i|23 



A l[(i,(A)|/,(A))(5(A)|,(A))|"'J 

By definition < FF < 1, and the effect of an imperfect matching of filters with signals translates into a reduction of 
the maximum SNR at which a source can be detected according to 

It has become standard in the gravitational wave data analysis literature to assume that a family of filters can be 
considered adequate if FF > 0.97, which corresponds to a decrease in detection rate by about 10%; the detection rate, 
in fact, scales as (FF)"^. Analogously, the match M is defined as 

„„.«,.,,.( "■""1"^)' \. (30, 

«-«..l|(A(A)|/i(A))(,(A)|5(A),l"=J 

where 6 and are the parameter vectors of the signal and template, respectively; A/o and A0o represent the 
mismatch in extrinsic parameters between signal and template. Notice that here we assume the use of the power 
spectrum as the detection statistic, which implies the need to search explicitly over the four angles 9]\[,(J)n,0l and 
(j)L. Instead of the power spectrum, one could actually consider the so-called F-statistic _21, 22], which requires the 
maximisation of the detection statistic using a discrete template bank only over the two angles 9^ and (f>N. The 
conclusions of this section are only marginally affected by our choice of detection statistic. 
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FIG. 3: The frequency dependence of the match. The plot shows the match A4, Eq. HHOfl . as a function of frequency, for four 
randomly selected sky positions and orbital plane orientations (dots: ^jv — 2.5, t^jv = 1.4, 9l = 1.2, and (f>L = 1-4; squares: 
On = 1.2, (l)N = 4.3, 6l = 0.36, and (J>l = 4.2; triangles: 9n = 3.1, (f)N = 0.81, Ol = 2.7, and <;/>£ = 0.79; circles: Gn = 1.8, 
(j)N = 3.6, &L — 1.9, and (/)_l = 3.6). The insert offers a zoom of the region around 10 mHz. The observation time is 1 yr. 



We aim to explore the FF for the family of filters constructed using the long wavelength approximation. The FF 
depends on the point in parameter space at which it is evaluated and its computation requires the maximisation over 
six parameters, cf Eq (|28|l . which in turn translates into a substantial computational burden. In order to investigate 
extensively the behaviour of templates based on the long wavelength approximation, while keeping the computational 
time at a manageable level, we have actually evaluated, over an extended portion of the whole parameter space, the 
match, Eq. ^U^, which represents a lower limit to the FF, and investigated the FF only for selected points. This 
strategy has the advantage of reducing substantially the computational burden, while still providing a reasonable 
approximation to the FF. We therefore set q = /il, cf Eq (|21|l . in Eq. (|30|l and H28|) and the signal is modelled 
according to Eq. (|15|) . All the results presented here assume that the observation time is T = 1 yr. 

The match and fitting-factor depend on both frequency and geometric parameters. In general we expect the fitting 
factor to be essentially unity for / <C /*, but to drop below 0.97 at some transition frequency whose exact value 
depends on ^^r, 0jv, Ol and <pL- In order to develop some intuition we begin by computing the match (we set 
& = 6 in Eq. (|30|l ) for four randomly chosen sky positions and orbital plane orientations as a function of frequency; 
the results are plotted in Figure |21 It is clear that the match falls below 0.97 only at about 10 mHz and is not a 
monotonic function of frequency; in fact there is no maximisation over the angular parameters and the impact of 
the transfer function depends on frequency in a non trivial way. Not surprisingly, the value of the match depends 
quite significantly on the position and orientation of a source in the sky. To investigate this dependence we select 
four fiducial frequencies (/q = 3 mHz, 5 mHz, 10 mHz and 30 mHz) and we compute the match for 1000 randomly 
chosen sky positions and wave polarisations at each of these frequencies |^. The results are summarised in Figure^ 
and Table Q] It is clear that in the frequency range of peak sensitivity of LISA (/ « 3 mHz) the long wavelength 
approximation is perfectly adequate for signal detection: all sources have a match (and a fortiori a fitting-factor) 
larger than 0.99. The frequency / = 10 mHz seems to mark the transition to the frequency range over which the exact 
LISA response is actually needed in order to recover the full SNR. Figure01and TableQlshow that at lOmHz the match 
for half of the sources is below 0.97, but always greater than « 0.91. In order to check whether the fitting-factor, 
the actual figure of merit we are interested in, shows a similar behaviour we have computed FF, Eq. (|28() . for the 
three sources that gave the three lowest values of the match, 0.914, 0.927 and 0.928. In all three cases, when the 
maximisation is carried out over the entire set of parameters, the fitting-factor is raised above 0.97, and yields 0.976, 
0.987 and 0.980, respectively. This is a very strong indication that FF is indeed > 0.97 at 10 mHz, and one can 
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FIG. 4: The distribution of the match. The plot shows the fraction of sources for which a given value of the match is attained: 
in each Monte Carlo simulation the match is computed for 1000 sources (all emitting at the same frequency) with random 
sky position and orientation of the orbital plane. Four values of the gravitational wave frequency are presented: /o = 3 mHz 
(dotted-dashed line), fo~5 mHz (solid line), /o — 10 mHz (dashed line) and /o = 30 mHz (dotted line). The observation time 
is T = 1 yr. The inserted panel shows the same results on a different scale for the x-axis in order to cover the whole range of 
values of the match for /o = 30 mHz. 

consider the long wavelength approximation still adequate at this frequency. However, at 30 mHz the picture changes 
dramatically and essentially all the sources yield a match which is below 0.97 - the largest match is 0.976 - with 
the lowest value being as small as 0.074. The mean and median are both « 0.5. In this case too we have computed 
the FF for a few sources. Placing the sources in order of ascending match we chose to find FF for the 1st, 2nd, 
300th and 700th entries, which gave matches of 0.074, 0.109, 0.432 and 0.558 respectively. As in the case of /o — 10 
niHz the maximisation over all the parameters increases quite substantially the actual value of the FF, in this case 
to 0.633, 0.614, 0.673 and 0.737 for the 1st, 2nd, 300th and 700th source, respectively. These results indicate that 
the threshold FF = 0.97, which we have set in order to consider the long wavelength approximation as suitable for 
detection, is not attainable for the large majority of sources. The outcome of our analysis is therefore clear: the long 
wavelength approximation, h ^ h^, Eq. I|21|l . can be safely used to construct detection templates up to /o — 10 mHz 
(for a fitting-factor > 0.97), after which it becomes progressively inaccurate, causing a severe reduction of detection 
rate; at / = 30 mHz the long wavelength approximation is by far inadequate. 



IV. PARAMETER ESTIMATION 



In this section we discuss the errors associated with the parameter measurements. The estimation of the minimum 
mean squared errors with which source parameters can be extracted in LISA observations of monochromatic sources 
has been investigated so far using the long wavelength approximation across the whole LISA observational window [7|, 
la, 13, llJI ■ The goal of this section is to explore to what extent the results reported in previous analyses are affected if 
one actually models the LISA output exactly. Results obtained using the long wavelength approximation are clearly 
correct at low frequencies (/ <SC /* « 10 mHz); however they become progressively less representative of what one can 
actually achieve during the mission in the high frequency portion of the LISA observational window. Here we establish 
the frequency at which the errors significantly differ depending on the approximation and quantify the discrepancy 
(as a function of frequency) with respect to previous work. We start by briefly recalling the general concepts and 
formulae regarding parameter estimation - we refer the reader to [3, Ha. Il9l |2J, |23 and references therein for more 
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/o 
(mHz) 


Match 
Min. Max. Median Mean Std. Dev. 


Fraction< 0.97 


3 


0.993 1.000 0.998 0.998 0.001 





5 


0.980 1.000 0.994 0.994 0.003 





10 


0.914 0.999 0.972 0.973 0.014 


0.437 


30 


0.074 0.976 0.506 0.515 0.183 


0.998 



TABLE I: Statistical summary of the values of the match obtained in each of the Monte Carlo simulations (same as Figure [IJ 
carried out over 1000 random source positions and orientations for a given gravitational wave frequency (shown in the first 
column of the table). The table shows the minimum (Min.) and maximum (Max.) value of the match, cf Eq. (I30L for the 
different values of the geometrical parameters 9n, (J^n, 9l, 4>l in each Monte Carlo simulation, and well as the (sample) median, 
mean and standard deviation (Std. Dev.) calculated over the 1000 values in the simulation. The last column shows the fraction 
of sources for which the match is < 0.97. 

details - and then present the resuhs of our analysis. 

In the hmit of large SNR, which applies to the vast majority of signals considered here, the errors AA associated 
with the parameters A that characterise h{t; A) follow a Gaussian probability distribution: 



p(AA) 



det(r) 
27r 



1/2 



g-iTj^AA^AA* 



(31) 



In Ea. H31() the matrix Tjk is known as the Fisher information matrix, which reads |3| 



M = 
^ jk - 



dh^'^ 



d\3 



a/i(') 



aA'= 



(32) 



the superscript i = 1, 1 1 labels the detector (as pointed out by Cutler in Q, using the three arms of the LISA 
instrument one can actually construct two observables, /i^^' and ft.*^^^', respectively, that correspond to two co-located 
interferometers one rotated by 7r/4 with respect to the other with uncorrelated noise.) The variance-covariance matrix 
is simply given by the inverse of the Fisher information matrix: 



E^*^ = (AA^' AA'') = 



r(') 



jk 



(33) 



The matrix S contains full information about the parameter errors and their correlations; in fact the diagonal elements 
of S represent the expected mean squared errors 



((AA^y) = E^'^' , 



(34) 



and its off-diagonal elements provide information about the correlations among different parameters through the 
correlation coefficients c-''^: 



c^' 



S-J* 



^^iJ^Ifc 



(-1 < d'' < +1) . 



(35) 



In the limit of high signal-to-noise ratio, E^-' pro vides a tight lower bound to the minimum mean-squared error 
((AA-')^), the so-called Cramer- Rao bound [l9ll25l |. Notice that the errors H34|l and the correlation coefficients H35|) 
depend on the actual value of the signal parameter vector A. For the case of observations with two or more detectors 



.W 



with uncorrelated noise, the Fisher information matrix is simply: Tjk = X^i, r^^. Two parameters we are interested 
in computing is the error associated with the position of a source in the sky (angular resolution) and the orientation 
of the orbital angular momentum. Following Q we define these by 



An 



N,L - 2n {(Acos6i^^i) (A^^^) - {A cos 6 nx A(l>Nxf] 



1/2 



(36) 



where N labels position and L labels orientation. The physical meaning of AnN,L is the following: the probability of 
N,L to lie outside an (appropriately shaped) error ellipse enclosing a solid angle Afl is simply e-^^/^^N,i,_ 
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b 


/ 






{^A/A)l/{^A/A) 






(mHz) 


Min. 


Max. 


Median 


Mean 


Std. Dev. 


Fraction < 1 


I 


3 


0.756 


1.642 


1.007 


1.049 


0.124 


0.40 


I+II 


3 


1.010 


1.052 


1.034 


1.033 


0.007 


0.0 


I 


5 


0.752 


2.483 


1.017 


1.127 


0.279 


0.358 


I+II 


5 


1.055 


1.127 


1.096 


1.094 


0.012 


0.0 


I 


10 


0.768 


6.535 


1.065 


1.511 


1.032 


0.303 


I+II 


10 


1.331 


1.541 


1.450 


1.444 


0.036 


0.0 


I 


30 


0.392 


44.164 


1.599 


4.949 


7.696 


0.219 


I+II 


30 


9.539 


36.435 


15.467 


16.769 


4.979 


0.0 



TABLE 11: Statistical summary of the values of {AA/ A) l / [^A/ A) obtained in each of the Monte Carlo simulations carried 
out over 1000 random source positions and orientations for a given emission frequency (shown in the second column of the 
table). The table shows the minimum (Min.) and maximum (Max.) value of (AA/A)_l/(AA/A) (see text for more details) 
for the different values of the geometrical parameters On, <Pn, Ol, 0l in each Monte Carlo simulation, as well as the (sample) 
median, mean and standard deviation (Std. Dev.) of {/S.A/A)l/{^A/A) calculated over the 1000 values in the simulation. 
The last column shows the fraction of sources for which {A.A/A)l/(/^A/A) < 1, that is the long wavelength approximation 
underestimates the errors in the measurement of the signal amplitude. The values are shown both for the case of observations 
carried out with a single interferometer (t — I) and with a pair of interferometers {i — I + II). The observation time is one 
year. Histograms for the distribution of {AA/A)l/{^A/A) are presented in Figure|3 



Each element of the Fisher information matrix, cf Eqs H26() and H32|l can be written as 

rj-fe^^ / 5,/iW(i)a,/.W(i)di, 

•JO J-oo 



(37) 



where dj = d/dX^ . Equivalently, the signal-to- noise ratio is given by 






(38) 



For the problem at hand, for which f ^ B ^dB/dt and / ^ B^^dBn/dt, cf Eq. (|^ and fTI)|l . one can actually 
simplify the full expression of the Fisher information matrix reducing it to 



1 f°" 
Tjk = ^ / [d,B{t)dkB{t) + B^t)d,x{t)dkX{t)] dt, 

•JO J-oo 

in the long- wavelength approximation - h — h^, ci Eq. (|21|l - and to 



(39) 



r.'- - i E 



So 



+ 



djBn{t)dkB,n + Bn{t)B„i{t)djXn{t)dkXin{t) COs(Xn - Xm) 



\Am{t)djAn{t)dkXm ~ An{t)djA,n{t)dkXn sin(x„ - Xm)dt 



(40) 



in the rigid adiabatic approximation, where h is given by Eq. H15|l . 

LISA parameter estimation for monochromatic signals strongly depends on the actual value of the signal parameters, 
in particular, emission frequency and location and orientation of a source with respect to the detector. This represents 
a large parameter space that one needs to explore in order to obtain meaningful results. We perform this exploration 
by means of Monte-Carlo simulations |23|. For each fiducial source we set /o and randomly select the geometrical 
parameters 9n, 4>n, Ol, (f'L (and the arbitrary initial phase 4>o) and compute the inverse of the Fisher information 
matrix for the cases where h is given by Eq. (|21|) and Eq. H15|l . The observation time is set to T = 1 yr. The signals 
are normalised in such a way that they produce the same signal-to-noise ratio in both cases. For each parameter 
A-', we can therefore evaluate the ratio (AA^)l/AA-', where the subscript 'L' indicates that the estimated error is 
computed using the long wavelength approximation for the LISA output; (AA-')l/AA'' is independent of SNR, and 
reflects simply the effect of structure of the transfer function on the estimation of the signal parameters. If the 
long wavelength approximation is indeed a good approximation for exploring the quality of LISA astronomy, then 
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FIG. 5: Distribution of the ratio between the errors associated with measurements of the signal amplitude using the long 
wavelength approximation and the rigid adiabatic approximation. The histograms show the fraction of sources (out of the 
1000, with random position and orientation, used in the Monte Carlo simulations) for which the amplitude can be measured 
with a fractional error characterised by a given value of the ratio (AA/yl)i/(Aj4/^), shown on the 2;-ajcis (see text for more 
details). Four histograms are shown, corresponding to different values of the gravitational wave frequency /o: 3 mHz (dotted- 
dashed line), 5 mHz (dotted line), 10 mHz (solid line) and 30 mHz (dashed line). One year of integration is assumed in the 
computation of the Fisher Information matrix. The top and bottom panels refer to measurements carried out with one and 
two detectors, respectively. A statistical summary of these results is presented in Table ITTI 



(AA-')l/AA^ = 1 . If (AA-')l/AA^ > 1(< 1) it means that results that have been presented in the literature so far 
overestimate (underestimate) the errors associated with parameter measurements and one can actually expect better 
(worse) quality astronomy with LISA. 

We have calculated the ratio (AA^)i/AA-' for four specific frequencies in the LISA observational window, /o = 
3 mHz, 5 mHz, 10 mHz and 30 mHz (the same values that we adopted in the study of match and fitting factor), for 
1000 randomly selected positions and orientations of sources in the sky. We have also considered both the case for 
one detector and for two detectors at 45° to each other. 

We present in detail results for three key parameters: amplitude, angular resolution of N and angular resolution of 
L in Figures 13 and [3 respectively, and Tables ^ IIIII and IIVI Each plot is a histogram showing the fraction of 
sources for which (AA-')l/AA^ (where j labels the relevant parameter) falls into each bin (for observations with one 
and two detectors). The tables provide a statistical summary of the outcome of each Monte Carlo simulation for the 
three chosen parameters. 

The qualitative behaviour of the results is very much consistent with what one would expect: as the frequency of 
the signal increases the effect of the LISA transfer function becomes more pronounced, and the errors on the source 
parameters are actually smaller than the ones obtained by using the simple long wavelength approximation, that 
is the ratios (A^/A)i/(Ayl/yl), (AilAr)L/Ai7Ar and (AriL)i/AJ7L become all greater than 1. This is clearly due 
to the structure of the signal recorded at the detector which is amplitude and phase modulated in a characteristic 
way that depends on the frequency and source position in the sky. In other words the actual LISA output provides 
more discriminating power than the one inferred from the simple long wavelength approximation. The values of 
the mean and median of (AA/^)i/(A^/A), (Arijv)L/AriAr and (Arii)i/Ar2i are all greater than 1 (i.e. the rigid 
adiabatic approximation provides errors smaller than the long wavelength approximation) and increase as the emission 
frequency increases (the higher the frequency the stronger the modulations, cf Figs. ^ and |2J|. Also the fraction of 
points for which the above ratios are larger than 1 increases. For / < 10 mHz the difference in results between the 
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FIG. 6; Distribution of the ratio between the errors associated with measurements of LISA's angular resolution using the long 
wavelength approximation and the rigid adiabatic approximation. The histograms show the fraction of sources that can be 
resolved in the sky with an angular resolution characterised by a given value of the ratio (Ar2]v)L/Af2jv, shown on the a;-axis 
(see text for more details). The parameters are the same as the ones in Figure |3 A statistical summary of these results is 
presented in Table UTTI 



L 


/ 

(mHz) 


(Afiiv)i/Af7iv 
Min. Max. Median Mean Std. Dev. Fraction < 1 


I 
I+II 


3 
3 


0.353 7.606 1.011 1.035 0.282 0.418 
0.560 10.666 1.036 1.062 0.366 0.185 


I 

I+II 


5 
5 


0.292 7.737 1.019 1.047 0.290 0.403 
0.561 11.832 1.099 1.127 0.408 0.380 


I 
I+II 


10 
10 


0.208 8.144 1.045 1.089 0.331 0.360 
0.676 17.225 1.446 1.490 0.602 0.010 


I 
I+II 


30 
30 


0.345 14.790 1.354 1.487 0.838 0.243 
6.808 259.199 18.908 20.586 12.545 0.0 



TABLE III: Statistical summary of the values of {AQn)l/^^n obtained in each of the Monte Carlo simulations. The 
parameters of the simulations are the same as the ones described in Table IITI We refer the reader to the former table for more 
details. Histograms for the distribution of (Anjv)L/Af2jv are presented in Figure|H| 



long wavelength approximation and the rigid adiabatic approximation is (on average) « 5 — 15%. For / > 10 mHz 
the discrepancy in the errors becomes much more pronounced. The difference in errors provided by the two different 
approximations is particularly strong when one considers the potential of LISA to behave as two interferometers with 
uncorrelated outputs: at 10 mHz and 30 niHz respectively the errors are actually 50% and a factor ^ 10 smaller than 
the ones obtained using the long wavelength approximation. However, it is important to stress - cf the spread in the 
histograms in Figures [SJ Eland [TJ and the minimum and maximum value of {AA/A)l/{AA/A), {AilN)L/AilN and 
(AJ7i)i/Ar2i in Table HTllIIII and |l VI - that if one considers any given source, even a.t fo — 3 mHz or 5 mHz right in 
the heart of the LISA sensitivity band where one would expect the long wavelength approximation to be completely 
adequate, the difference in the value of the errors can range from « 30% to a factor « 10. Moreover, depending on 
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FIG. 7: Distribution of the ratio between the errors associated with the measurements of the orientation of the binary orbital 
plane using the long wavelength approximation and the rigid adiabatic approximation. The histograms show the fraction of 
sources whose orbital plane orientation, described by the unit vector L (cf the Appendix) can be resolved within a solid angle 
characterised by a given value of the ratio {AQl)l/AQl, shown on the x-axis (cf text for more details). The parameters are 
the same as the ones in Figures |^ and |S| A statistical summary of these results is presented in Table ITVI 



t 


/ 

(mHz) 


(Ar!i)i/Ar!i 

Min. Max. Median Mean Std. Dov. Fraction < 1 


I 
I + II 


3 
3 


0.422 1.384 1.017 1.043 0.077 0.295 
0.995 1.086 1.033 1.032 0.010 0.006 


I 

I + II 


5 
5 


0.915 2.312 1.036 1.110 0.192 0.260 
1.011 1.154 1.094 1.092 0.016 0.0 


I 
I + II 


10 
10 


0.825 5.187 1.105 1.440 0.777 0.235 
1.062 1.614 1.442 1.435 0.053 0.0 


I 
I + II 


30 
30 


0.247 61.811 1.856 4.746 7.429 0.211 
5.988 52.498 15.150 16.941 6.567 0.0 



TABLE IV: Statistical summary of the values of {AflL)L/A^L obtained in each of the Monte Carlo simulations. The 
parameters of the simulations are the same as the ones described in Tables ITTl and 11111 We refer the reader to the former tables 
for more details. Histograms for the distribution of {AQl)l/AQl are presented in Figure|7| 



the source geometrical parameters, the actual errors (the ones derived using the rigid adiabatic approximation) could 
be larger than the ones predicted using the long wavelength approximation. For detailed studies of LISA astronomy 
it is therefore important to consider the real response of the instrument even in the mHz range. 

We would like to conclude this discussion by pointing out another subtlety. So far we have concentrated on the 
comparison of the lower bounds on the minimum mean squared errors estimated using two different expressions of 
the detector output. The conclusion is that the full transfer function needs to be included in order to obtain accurate 
predictions of the quality of LISA astronomy. If one sticks to the long wavelength approximation as a model of the 
detector output, such predictions become increasingly less accurate as the frequency of the putative source increases. 
We have therefore argued that the errors are (on average) actually smaller than the ones reported so far in the 
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literature. However, the actual errors with which LISA will be able measure source parameters do not necessarily 
decrease at higher frequencies: in fact, for a given distance, position and orientation in the sky of a monochromatic 
source, when the frequency increases also the LISA noise spectral density becomes larger (which in turn degrades the 
quality of measurements). There is therefore a competition of effects, and it is simple to show that eventually the 
signal-to-noise degradation "wins" over richness of signal structure. 

V. CONCLUSIONS 

We have explored the effect of the LISA transfer function on signal detection and parameter estimation for monochro- 
matic sources. We have shown that for signal detection, although the long wavelength approximation is not a faithful 
representation of the Michelson interferometer output for f ^ 3 mHz, it can actually be safely used in constructing 
templates for monochromatic signals up to 10 mHz. The fitting-factor, in fact, always exceeds 0.97 in this frequency 
regime. For parameter estimation, on the other hand, the effects of the frequency and position dependent transfer 
function already become significant around 3-to-5 mHz and the use of the long wavelength approximation is totally 
inappropriate for / ^ 10 mHz: the estimate of the measurement errors can be quite different, ranging from 5% to a 
factor ~ 10 (or more), from the ones obtained using the low frequency approximation depending on the source param- 
eters and the observational mode (LISA as a 2 or 3 arm instrument). Our analysis is limited in two main respects. In 
the study of signal detection we have explored the parameter space quite extensively for the calculation of the match, 
but not for the fitting factor; this is entirely due to the computational constraints, but can be alleviated in the future 
by accommodating longer runs. Secondly, both for parameter estimation and signal detection, we have restricted 
attention only to monochromatic signals. If the source chirp mass is, say, 1 Mq and the observation time 1 yr, the 
radiation becomes linearly chirping at / « 5 mHz. This introduces an additional parameter (the first time derivative 
of the frequency) in the signal gravitational wave phase. The transition to a linearly chirping regime depends on the 
source chirp mass and is at higher frequencies for white dwarf binary systems. Therefore, depending on the source 
mass that is considered, the analysis presented in this paper might require the addition of one parameter in the signal 
waveform. However, it is clear that the results presented here would not be affected in any significant way, and the 
general conclusions regarding the effect of the transfer function would apply. We plan to return to these issues in a 
future paper. 
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APPENDIX A: FORMULAE 

Here we present explicit analytical expressions for the rigid adiabatic approximation h(t), see Eqs. (I15ll - (|18|) . and the 
long wavelength approximation hL{t), see Eqs. (I^ - (|^ . to the LISA detector output that are discussed in SectionllTl 

It is convenient to introduce two Cartesian reference frames (cf, e.g., \^)'. a "barycentric" frame {x, y, z) tied to the 
Ecliptic and centred in the Solar System Barycentre, with z perpendicular to the Ecliptic, and the plane (a;, y) in the 
Ecliptic itself; a detector reference frame (a;', y', z'), centred in the LISA centre of mass and attached to the detector, 
with z' perpendicular to the plane defined by the three arms and the x' and y' axis defined so that the unit vectors 
\j {j = 1,2,3) along each arm read 



L- = cos 



12 + 3(^'-^^ 



T2 + 3(^'-^^ 



- sin aj{t) cos<I>(i) — cos a j{t) sin<I>(t) 



- smaj{t) sin<I>(t) + cos a j{t) cos$(f) 



— smaj(i) 



(Al) 



where aj(t) increases linearly with time, according to 



aj{t) = n^t - (j - l)7r/3 + ao 



(A2) 
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In the previous expression ao is a constant specifying the orientation of the arms at the arbitrary reference time t = 0, 
and Uq = 27r/lyr. In the Echptic frame the motion of LISA's centre-of-mass is described by the polar angles 



e^f 



and the normal to the detector plane z' precesses around z according to 



V3 



cos $(i)x + sin $(t)y] 



(A3) 



(A4) 



We shall follow the convention that primed and unprimed quantities refer to the frame attached to LISA and the Solar 
System Barycentre, respectively. The geometry of a binary system with respect to LISA is described by the two unit 
vectors N and L, where the former identifies the source position in the sky and the latter defines the direction of the 
orbital angular momentum (spins are negligible for sub-solar mass binary systems and L can be therefore regarded as 
constant). With respect to an observer on LISA, i.e. the reference frame {x' , y' , z'), the polar coordinates of N and L 
are {d'j^ , (p'pf) and (6*^,0^), respectively. We also define ip' as the time dependent polarisation angle. The angles 0^, 



i/jY and ip' can be written as a function of the angles measured with respect to the solar system barycentre as 

1 



tan ■(/'at = 



■ cos t/]\j 

■K 

'^+12 



— sin 6m cos($(t) - 07v) , 



tan 



VS cos Om + sin 9m sin(<i>(i) — (J)m) 



2sin0Arsin($(i) - 0Af) 



L • z' - (L • N) (z' • N) 
N • (L X z') 



where 



Ej = n(^t ■ 



n 
12 



rO- 



1) + Sn 



and Sq sets the orientation of 1^ at t = 0. 

The relevant scalar and vector products entering the definition of ip' are 



L • z' = - cos 9r 
2 



— sin6lLCOs($(t)-0L), 



N 



L-N 
Lx z' 



cos 9l cos 9n + sin 9l sin 9n cos(0l — ^jv) , 
- sin 9]y sin 9^ sin((/)i — (f)^) 



COS $ (t) ( COS 9l sin 9n sin (p^ 

sm $ 

2 



COS 9n sm Ul sm (p^ 



W(' 



COS Un sm c/i cos cpL — cos 0^ sm 0^ cos 0jv 



The scalar products Ij • N (j = 1,2,3) entering the expression of the LISA transfer function, Eqs (|S 
and Hlla|) - (|lld|) can be computed in a straightforward way using Eqs. (|A1I) . lASp and (|A7ll . 

The wave polarisation amplitudes, that enter the expressions H16|l . (|18|1 . (|22|l and H24() . are given by 



A+{t) 



Ax(t) = -4 



D 
D 



1 



(l.n)^ 



L-N {-Kf) 



irrf) 



2/3 



2/3 



(A5) 
(A6) 

(A7) 
(A8) 



(A9) 
(AlO) 



(All) 
(tTTTali-inMli 



(A12a) 
(A12b) 



where A4 is the source chirp mass, defined as A^ = m ' /i ' - m and /i are the total and reduced mass of a binary 
systems, respectively -, / is the GW frequency and D is the (luminosity) distance. The GW phase in Eqs. H17|) 
and l|^ is 0gw(*) = 27r/ot — (po- 
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The Doppler phase modulation, in Eqs. (|14() . H17() and H23|l . is 

ifuit) = 2nR(s f sinON cos($(i) - (I)n) , (A13) 

where Rm — 1 AU. 
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